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The electronic properties, exchange interactions, finite-temperature magnetism, and transport 
properties of random quaternary Heusler Ni2MnSn alloys doped with Cu- and Pd-atoms are studied 
theoretically by means of ab initio calculations over the entire range of dopant concentrations. While 
the magnetic moments are only weakly dependent on the alloy composition, the Curie temperatures 
exhibit strongly non-linear behavior with respect to Cu-doping in contrast with an almost linear 
concentration dependence in the case of Pd-doping. The present parameter-free theory agrees quali- 
tatively and also reasonably well quantitatively with the available experimental results. An analysis 
of exchange interactions is provided for a deeper understanding of the problem. The dopant atoms 
perturb electronic structure close to the Fermi energy only weakly and the residual resistivity thus 
obeys a simple Nordheim rule. The dominating contribution to the temperature-dependent resis- 
tivity is due to thermodynamical fluctuations originating from the spin-disorder, which, according 
to our calculations, can be described successfully via the disordered local moments model. Results 
based on this model agree fairly well with the measured values of spin-disorder induced resistivity. 

PACS numbers: 71.23.-k,72.25.Ba,75.10.Hk,75.30.Et 



G 

O 

o 



> 
o 

O 

od 
O 

o 



X 



I. INTRODUCTION 

Heusler alloys were first studied by the German 
chemist Friedrich Heusler in 1903, starting with the or- 
dered alloy Cu2MnSn. Because of their interesting phys- 
ical properties they have been studied intensively in the 
past as well as more recently^ The most widely stud- 
ied Heusler alloys are those with the formula Ni2MnZ 
(Z=Sn,Ga,In) and Co 2 XY (X=Mn, Fe; Y=A1, Si). The 
former group is of interest because of potential technolog- 
ical applications based on their magnetic shape memory 
effect. 2 the magnetocaloric effect,' 5 and the recently ob- 
served giant (negative) magnetocaloric effect X The latter 
group holds the promise of application in spintronic de- 
vices, thanks to their halfmetallicity at room temperature 
and above, lattice constant matching with the III-V semi- 
conductors, and large bandgaps. Large tunneling mag- 
netoresistance was measured recently in Co2MnSi/Al- 
0/Co2MnSi magnetic tunneling junctions, 5 

Structurally, most Heusler alloys crystallize in two 
different cubic phases, having either the L2i (X2YZ) 
or the Cib (XYZ) symmetry. They can be best visu- 
alized as being composed of four interpenetrating fcc- 
sublattices, shifted along the body diagonal in the or- 
der X-Y-X-Z or X-Y-E-Z, where E denotes the empty, 
i.e. unoccupied, sublattice. 1 An important feature of 
Heusler alloys is the presence of chemical or substi- 
tutional disorder. It is often the non-stoichiometric 
composition with respect to the ideal systems such 



as Ni2MnSn or N^MnSb which interpolates between 
the L2i Heusler and the Cib semi-Heusler alloys. 6 Ex- 
amples are the magnetic shape memory alloys of the 
type Ni2Mni-|_ a; Siii_ a; (Mn-nonstoichiometry), or the 
Ni2- x MnSb (Ni-nonstoichiometry) alloys. 

In addition to the chemical disorder due to nonstoi- 
chiometry, a native chemical disorder exists even in 'ideal' 
ordered alloys X2YZ and XYZ. In particular, halfmetallic 
Heusler alloys are very susceptible to such native disor- 
der, a typical example being the Co2MnSi alloy which 
exists in the B2-like structure due to Mn-Si disorder. 

Finally, there are complex quaternary alloys like the 
semi-Heusler (Ni,Cu)MnSb alloys^ and closely related 
Heusler alloys like (Ni,Cu) 2 MnSn 9 or (Ni,Pd) 2 MnSn^ 
with disorder on the X-sublattice. The magnetic prop- 
erties of these alloys are simplified by the fact that the 
Mn-sublattice is effectively nonrandom (although there 
can be some magnetic disorder in a certain concentra- 
tion range) i 7 - In addition, Mn-atoms carry essentially all 
of the system magnetic moment. However, despite these 
simplifying features the magnetic, thermodynamic, and 
transport properties of these alloys are not easy to under- 
stand. Several factors: disorder on sublattices neighbor- 
ing the unperturbed Mn-sublattice, varying carrier con- 
centration, the presence of atoms with different degrees 
of d-electron localization and thus different levels of hy- 
bridization with Mn-atoms, all contribute to the com- 
plexity of the problem. A reasonably successful interpre- 
tation of the results of experiments in Refs.[9lfl0j was pro- 
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vided by Stearn a 11 ' 12 using a qualitative model. Accord- 
ing to this, magnetic behavior of Heusler alloys is con- 
trolled by three types of magnetic interactions, namely: 
(i) interaction between extended s-like electrons and lo- 
calized d-electrons mediated by the sd-hybridization, (ii) 
interaction between localized and delocalized (itinerant) 
d-electrons, and (iii) the superexchange mediated by sp- 
element on the Z-sublattice (Sn in the present case). The 
dominating interaction is that between the localized and 
itinerant d-electrons and as such the value of the Curie 
temperature T c of Heusler alloys is essentially controlled 
by the amount of itinerant d-electrons. This, in turn, is 
related to the localization of the corresponding d-orbitals 
of the X-atoms. All of the above mentioned interactions 
are naturally present in the first-principles description of 
magnetism in the framework of the spin density func- 
tional theory (DFT). In general, it is very difficult to 
separate out the individual contributions from the calcu- 
lated total interactions. 

In this work we employ state-of-the-art electronic 
structure calculations to understand the properties of 
disordered quaternary Heusler alloys (Ni,Cu)2MnSn and 
(Ni,Pd)2MnSn. The exchange interactions are computed 
from the ab initio electronic structure, and then used to 
estimate the T c using various approximations. Such a 
program has been carried out in recent years for a num- 
ber of conventional (ordered) Heusler alloys. Notable are 
the pioneering work of Kiibler- 3 and extensive studies of 
electronic and magnetic properties of Heusler alloys, in- 
cluding estimates of the Curie temperature, by Kiibler, 14 
Galanakis, 15 Picozzi^ and Sasioglu and Sandratskii. 17 
The study of disordered Heusler alloys is, however, much 
more involved from a theoretical standpoint because of 
the presence of disorder violating the translational sym- 
metry. Although in some cases, such as the Heusler al- 
loys with sp-disorder on the Z-sublattice, it is possible 
to use conventional band structure methods (the super- 
cell approach),— more sophisticated methods, typically 
using the coherent potential approximation (CPA) im- 
plemented within DFT formalism (DF-CPA), need to be 
cmploye d 19 i 20 to treat disorder in general. This is partic- 
ularly true for systems with general compositions such as 
the magnetic shape memory alloys or the Heusler alloys 
in narrow concentration ranges where abrupt changes 
in physical properties are known to occur (e.g. around 
the austenite-martensite transition). So far, such stud- 
ies have been typically limited to the electronic structure 
and simple magnetic properties, e.g. magnetic moments. 
The first attempts to study thermodynamical properties 
(including T c ) of random alloys have appeared^ only 
recently. Particularly worth mentioning is the extensive 
study of the semi-Heusler alloys (Cu,Ni)MnSb 7 compris- 
ing the electronic structure, magnetic, thermodynamical, 
and transport properties of random alloys in the frame- 
work of a unified DFT description*! 

The aim of the present work is to study electronic, 
magnetic, thermal, and transport properties of two re- 
lated disordered Heusler quaternary systems based on the 



reference alloy ISfisMnSn, namely the (Ni,Cu)2MnSn and 
the (Ni,Pd)2MnSn systems. Both systems have magnetic 
moments around 4 /^b over the whole concentration range 
but dramatically different concentration dependence of 
their T c : strongly nonlinear in the former case and al- 
most linear in the latter. The T c of end-point alloys, 
namely X2MnSn, with X=Cu, Ni, and Pd, have decreas- 
ing values in this order. Thus they seem to obey the 
criterion advanced by Stearns^ for the dependence of T c 
on the localization of d-electrons of the element X. Lo- 
calization of the d-electrons increases left to right across 
a transition metal row, and decreases across a column as 
we go from 3- to 4- and then to 5-d. Thus Cu d-electrons 
are more localized than Ni d-electrons, which are more 
localized than Pd d-electrons, explaining the differences 
in T c of the three systems in the order mentioned above. 

Transport properties of Heusler alloys have not been 
studied extensively. Some measurements of the tempera- 
ture dependence of resistivity have been carried outj 21 i 22 
To our knowledge, no theoretical studies of transport in 
Heusler alloys have appeared so far. We attempt to fill 
this gap, using a simplified approach. 



II. FORMALISM 

The electronic structure calculations were performed 
using the tight-binding linear muffin-tin orbital (TB- 
LMTO) scheme 23 in the framework of the local den- 
sity approximation (LDA). The effect of substitutional 
disorder on the X-sublattice (either Ni-Cu or Ni-Pd) is 
described by CPA formulated in the framework of the 
TB-LMTO Green's function method^ The same atomic 
sphere radius was used for all the constituent atoms, and 
lattice constants were taken from experiments^^ The 
calculations employed an s,p, d, /-basis. For the parame- 
terization of the local density functional the Vosko-Wilk- 
Nusair exchange-correlation potential 2 ^ was used. 

The thermodynamical properties of the system are as- 
sumed to be given by a classical Heisenberg Hamiltonian, 

H c ff — — Jij e, • ej , (1) 

id 

where i,j are site indices, is the unit vector pointing 
along the direction of the local magnetic moment at site 
i, and Jij is the exchange integral between sites i and 
j. The exchange integrals, by construction, contain the 
atom magnetic moments, their positive (negative) val- 
ues being indicative of ferromagnetic (antiferromagnetic) 
coupling. 

We evaluate exchange integrals in Eq. (flj using a two- 
step model ) 26 i 27 where the band energy is equated to the 
Heisenberg form and then expressed via multiple scat- 
tering formalism based on the TB-LMTO-ASA Green's 
function in terms of the moments directed along and 
ej. The reference state for this calculation is chosen 
to be the disordered local moment (DLM) stated Such 
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a choice was recently suggested for the study of semi- 
Heusler alloys^ as it has some advantages over the con- 
ventional choice of the ferromagnetic reference state.— 
The DLM state is closer to the state at which the mag- 
netic transition occurs, compared with the state with a 
global magnetization. There is no preferred magnetic 
configuration assumed, and there are no induced mo- 
ments for the DLM state (see recent discussion of the 
problem of induced moments in the Heisenberg model 
in Ref|29j). The DLM reference state was successfully 
used recently in the study of magnetic overlayers on non- 
magnetic substrates^ Ideally one would hope the result- 
ing magnetic properties to be robust, i.e. (almost) inde- 
pendent of the assumed reference state. In practice, some 
dependence on the reference state is unavoidable, and in 
some cases the DLM reference state has been found to 
provide better estimates of the Curie temperatures. 

We determine the Curie temperature corresponding to 
the effective Heisenberg model in Eq. (TT|) by making use 
of the random-phase approximation (RPA). For compar- 
ison, we also include results obtained in the framework 
of the the mean-field approximation (MFA). The RPA- 
Curie temperatures are known to be close to those ob- 
tained from Monte-Carlo simulations <2! 

We have used interactions up to ~ 4 lattice constants 
which arc found to be sufficient to achieve reasonably 
converged results. While it is not a problem to include 
even more distant interactions in the evaluation of the 
RPA-Curie temperature, problems arise in connection 
with integration of the inverse lattice Fourier transform 
of real-space exchange integrals over the Brillouin zone 
which becomes a rapidly oscillating function for distant 
shells. As a result, the RPA Curie temperature oscillates 
as a function of the shell number included in the above 
mentioned lattice Fourier transform. It should be noted 
that such problems do not arise for the MFA or while 
dealing with semi-Heusler alloys^ where the exchange 
integrals are strongly damped in real-space due to their 
halfmetallicityi^ We have adopted the approach proposed 
m RefiH and used in the context of Heusler alloys in 
Ref@. It relies on using a set of exponential damping 
parameters to compute the T c , thereby reducing the os- 
cillations and finally extrapolating the results to the zero 
damping case. This approach was successfully applied to 
the evaluation of the stiffness constant of ferromagnets 
in real-space, where the problems associated with using 
only a limited number of shells are even more severed 
The above-mentioned oscillations are even stronger for 
the ferromagnetic reference state. 

Without an external magnetic field, there are essen- 
tially three different contributions to the resistivity of 
magnetic alloys: (i) phonon scattering, (ii) magnetic scat- 
tering due to thermodynamical fluctuations, which are 
largest close to the Curie temperature^ and (iii) the 
residual resistivity due to the presence of chemical disor- 
der on the Ni-Cu and Ni-Pd sublattices. 

In this study we limit ourselves to the residual resistiv- 
ity, and also present a simplified treatment of the resistiv- 



ity due to magnetic (spin-disorder) scattering. The resid- 
ual resistivity is determined by the linear-response theory 
as formulated in the framework of the TB-LMTO-CPA 
approach using the Kubo-Greenwood formula^ i.e., on 
the same formal footing as used for the determination 
of the exchange integrals. This approach, formulated 
for the multi-sublattice case, allows us to include both 
the substitutional disorder on the Ni-sublattice as well 
as magnetic disorder on the Mn-sublattice on an equal 
footing^ The disorder-induced vertex corrections are in- 
cluded in the formalism^ The spin disorder is described 
here in the framework of the DLM, and treated formally 
as 'substitutional' disorder via CPA. We refer the reader 
to Ref|H for details on the implementation of the DLM 
in the framework of the TB-LMTO formalism. 



III. RESULTS AND DISCUSSION 

In this section we present results for the electronic, 
magnetic, and transport properties of (Ni,T)2MnSn 
(T=Cu, Pd) alloys over a broad range of concentra- 
tions and compare our results with available experimen- 
tal data. 



A. Density of states 

Fig. [T] shows the calculated local densities of states 
(LDOSs) of nonmagnetic X 2 MnSn (X=Cu, Ni, Pd) al- 
loys. Several conclusions follow immediately: (i) all 
DOSs have a pronounced peak at the Fermi energy (£"f), 
indicating instability against formation of the ferromag- 
netic state (Stoner criterion), (ii) The dominant contri- 
bution to the peaks is from the Mn-LDOS, which sug- 
gests that the magnetism is primarily due to the Mn- 
sublattice. (iii) The height of Mn-LDOS(£ F ) (which is 
~ total DOS(-Ef)) is approximately the same in all cases, 
indicating similar magnetic moments for all these com- 
pounds (see Fig. [3] below). (iv)The increasing bandwidth 
of the d-states dominating local X-LDOS in the order 
Cu-Ni-Pd is clearly seen. This is indicative of decreasing 
localization (increasing itinerancy) of the d-electrons, re- 
sulting in decreasing T c in the same order, according to 
the model proposed by Stearns*^ (v) The Cu- and Ni- 
bands are well separated in energy resulting in a strong 
diagonal disorder, while their bandwidths are approxi- 
mately the same (weak off-diagonal disorder). On the 
contrary, Ni- and Pd-bands differ mainly in their widths 
(off-diagonal disorder). In all cases, however, X-bands 
are well separated energetically from the Fermi energy so 
that disorder on the X-sublattice influences states at E-p 
only weakly. Hence, the residual resistivity (due to chem- 
ical disorder) should be in the weak-scattering regime. 

The spin-polarized LDOS for Ni 2 MnSn, 
(Ni5o,Cu5o)2MnSn, and Cu2MnSn are shown in 
Fig. [5] (top panel). The following conclusions can be 
drawn: (i) Spin-polarized Cu-LDOSs in Cu2MnSn are 
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FIG. 1: Total and component resolved densities of states, per formula unit and spin, for nonmagnetic Heusler alloys: (a) ordered 
Cu 2 MnSn, (b) ordered Ni 2 MnSn, and (c) ordered Pd 2 MnSn. 



almost identical, indicating practically no polarization, 
while some small polarization is seen for corresponding 
Ni-LDOS in N^MnSn. Similar conclusions are also valid 
for Cu- and Ni-LDOSs for the equiconcentration case 
(Fig. [5Jd). (ii) The total LDOS is smoothed by strong 
level disorder in the equiconcentration alloy, (iii) Due 
to large level splitting (large local magnetic moment) 
the minority Mn-bands in all cases hybridize very little 
with the X-bands. On the contrary, such hybridization, 
compared to the nonmagnetic case, is strong for the 
majority bands even for Cu2MnSn. (iv) The majority 
and minority states at the Fermi energy behave differ- 
ently (corresponding DOSs have different curvature). 
This results in different Fermi surface geometry for 
these bands (see e.g. Refl37T). The results for Ni2MnSn, 
(Ni 5 o,Pd5o) 2 MnSn, and Pd 2 MnSn are shown in Fig. [2] 
(bottom panel). There are some differences in this 
case: (i) The carrier concentration is the same through 
the entire alloy system. As a result, the position of 
the minority Mn-band with respect to the Fermi level 
stays fixed. This is different from the (Ni,Cu)2MnSn 
alloy system where the carrier concentration increases 
with Cu-content. (ii) We observe strong hybridization 
of Ni- and Pd-bands on the X-sublattice, which differ 
mostly in their widths, (iii) One thus expects weaker 
site off-diagonal disorder and, consequently, also lower 
residual resistivity (see Section UlI Dl below) . 



B. Magnetic moments 



Experimental and calculated average moments per for- 
mula unit (f.u.) are presented in Fig. [3J together with 
calculated local moments on Mn-sublattices. There is an 
overall good agreement between calculated and measured 
average moments which are essentially concentration- 
independent and have values around 4 /j,b- For the stoi- 
chiometric alloys X2MnSn (X = Ni, Cu, Pd), very simi- 
lar values were obtained in a previous study by the TB- 
LMTO-ASA method^ Mn-sites carry almost all the mo- 
ment per f.u. and this is particularly true for Cu-rich 
alloys. Although the calculated magnetic moments agree 
with experimental values within a few percent, even bet- 
ter agreement can be obtained (e.g. for Ni2MnSn) when 
calculations are performed using the generalized gradient 
approximation (GGA)i£ instead of LDA. The weak con- 
centration dependence of the averaged magnetic moment 
can be understood from Fig. [2] The large intra-atomic 
exchange splitting of Mn atoms and the small hybridiza- 
tion of minority Mn-d orbitals with the X-atom orbitals 
lead to the full occupation of majority Mn-d orbitals, to 
negligible spin polarization of X- and Sn-atom orbitals, 
and to composition-independent occupation of minority 
Mn-d orbitals. These features result in the nearly con- 
stant Mn-moment and the total alloy magnetization. The 
local moments on Ni-sites have values in the range (0.14, 
0.19) hb for both systems and all concentrations. The 
local Cu moments are very small (0.01—0.02 /zb), while 
local Pd-moments are almost concentration-independent, 
with a value ~ 0.1 /Ub- Finally, local moments on the Sn- 
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FIG. 2: Total and component resolved densities of states, per formula unit and spin, for ferromagnetic Heusler alloys. Top panel: 
(a) ordered Ni2MnSn, (b) disordered (Niso,Cu5o)2MnSn, and (c) ordered Cu2MnSn. Bottom panel: (a) ordered Ni2MnSn, (b) 
disordered (NiscPdso^MnSn, and (c) ordered Pd2MnSn. Majority (minority) densities of states are shown in upper (lower) 
parts of each figure. 



sublattice are small and negative, in the range (-0.02, 
-0.04) hb. 



C. Exchange interactions 

1. Calculated results 

The exchange integrals are one of the important char- 
acteristics of the magnetically polarized state. We have 
determined them as a function of the alloy composition 
using the DLM reference state. It should be noted that 



for the DLM reference state the only non-zero exchange 
interactions are those between the Mn sites. The Mn- 
moments are very rigid and their values in the ferromag- 
netic and DLM states almost coincide with each other. 
The leading exchange integrals in (Nii_ a ,Cu :c )2MnSn 
quaternary Heusler alloys are plotted as a function of 
the Cu-concentration in Fig. HJi. The Mn-sublattice is 
only indirectly perturbed by the alloy disorder on the X- 
sublattices, via the hybridization of Ni- and Cu-atoms 
with the Mn-sublattice. Such a hybridization is weak 
for Cu-atoms, but strong for the Ni-atoms (see Figs. Q] 
and [2J . In addition to Mn-X hybridization, which varies 
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FIG. 3: Averaged magnetic moment (full circle) and local 
magnetic moments of Mn-atoms (diamonds) as a function 
of composition for (Ni,Cu)2MnSn alloys. Related results for 
(Ni,Pd)2MnSn are shown as inset. Experimental data^^ are 



denoted by empty circles. 



with the alloy composition, there is also an increase of 
the total electron concentration with the Cu-content due 
to different valency of Ni- and Cu-atoms. The resulting 
concentration dependence is quite complex, with the first 
three exchange integrals being dominant. While the first 
two are ferromagnetic (FM) and increase with the Cu- 
content, the third one decreases with Cu-concentration 
almost linearly and changes its sign from FM to antifer- 
romagnetic (AFM) at about 20% of Cu. It should be 
noted, however, that the 4th and in particular the 6th 
neighbor interactions are also non-negligible. To better 
understand the contribution of the various shells to the 
magnetic properties we have multiplied the exchange in- 
tegrals by their degeneracies (N), i.e., by the number of 
the equivalent atoms in the corresponding shells. This 
is presented in Fig. Such terms provide the actual 
contribution of a given shell to the MFA expression of 
the Curie temperature! 26 ' 27 With the multiplication by 
the degeneracy factors, the contribution of the first shell 
remains the largest one, and the contributions of the 2nd 
and 6th shells are suppressed compared to that of the 
3rd shell. We have also shown the contribution of the 
7th shell, which is non-negligible due to its large degen- 
eracy, even though individual atom contribution is very 
small. 

The above results confirm the assumptions of a the- 
oretical model of Stearns^ which asserts that at least 
three independent magnetic interactions are needed in 
order to explain the magnetic properties of Heusler al- 
loys with late 3 d and 4 d elements on the X-sublattices. 
This is a probable cause for the failure of a first-principles 
study 13 in which the authors tried to estimate the Curie 
temperature from just two exchange integrals, based on 
the energy differences between the FM state and two dif- 



ferent AFM states (AF I and AF II with antiparallel 
alignments of spins along the [fOO] and [fff] directions). 
A bad estimate of the Curie temperature for such a model 
for Cu2MnSn can be attributed particularly to the large 
contribution of the 3rd shell (see Fig. |4Jd), which was ig- 
nored. 

In Fig. [5] we present the results of a similar study for 
(Nii-^Pd^MnSn Heusler alloy. There are some impor- 
tant differences with respect to the (Ni,Cu)MnSn alloy: 
the number of valence electrons remains unchanged upon 
alloying, Ni and Pd being neighbors in the same column 
of the Periodic Table. As already mentioned, the effect of 
disorder is smaller, as it has predominantly off-diagonal 
character. The leading interactions are FM, but some 
interactions, e.g. those from the 4th and 6th shells, are 
AFM (see Fig. [3Jl) . Upon multiplying the exchange in- 
tegrals by their shell degeneracy (Fig. [SJd), the contribu- 
tions from the 1st and 3rd shell are seen to be dominant, 
while the contributions from the 2nd, 5th, and 7th shells 
(FM) and those of 4th and 6th shells (AFM) tend to com- 
pensate each other in the MFA sum. The contributions of 
all shells are either almost concentration-independent or 
decrease slightly with the Pd-content, resulting in a linear 
decrease of the Curie temperature with Pd-concentration 
(see below). 

There is a natural question as to how many shells are 
needed for a reliable estimate of the Curie temperature. 
We have just seen that the smallness of exchange inte- 
gral is not a sufficient condition, as it can be offset by 
a large degeneracy of the shell. We have insured that 
the weighted sum of exchange integrals as a function 
of the shell number achieves its saturation value. For 
halfmetallic alloys with exponentially damped exchange 
integrals^— this condition is fulfilled for a relatively small 
number of shells. In transition metal alloys with strong 
level disorder (see e.g. RefjH) the situation is similar, 
because the disorder influences the states at the Fermi 
energy, which are the relevant states dictating the spa- 
tial dependence of exchange integrals^ However, for the 
Heusler alloys considered here, the effect of disorder at 
the Fermi energy is not as strong as in halfmetallic alloys. 
Hence a larger number of shells needs to be included. We 
illustrate the situation in Fig. [6j where we plot the mean- 
field Curie temperature for three typical compositions in 
(Nii-^Cu^MnSn alloy. We note that the mean-field 
value of the Curie temperature is directly proportional 
to the sum of the exchange interactions We observe 
that for the shell distance d < 4 a, where a is the lattice 
constant, the results are reasonably saturated. We have 
thus used all neighbors up to 4 a and applied the extrapo- 
lation procedure briefly described in the Introduction. 31 
This procedure is necessary, as the RPA T c values os- 
cillate more strongly than the MFA values due to the 
contribution from the neighborhood of the zone center, 
where the contribution of the distant shells plays an im- 
portant role. The estimated error of this procedure is a 
few percent of the absolute value of the Curie tempera- 
ture. 
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FIG. 4: Concentration dependence of exchange integrals in (Nii-^Cu^MnSn Heusler alloys as estimated from the disordered 
local magnetic reference state: (a) the first six leading (Mn,Mn)-exchange integrals, and (b) the same but with exchange 
integrals multiplied by their shell degeneracies. We have also added the 7th shell term with a large degeneracy. 
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FIG. 5: The same as in Fig. g] but for (Nii-^Pd^MnSn Heusler alloys. 



Calculated results for (Ni,Cu)2MnSn are presented and 
compared with experiments^ in Fig. [Jj 

The system (Nii—^Cua^MnSn exhibits different be- 
haviors in two concentration regions: for x € (0.0, 0.3) 
the Curie temperature stays almost constant with a small 
minimum at x w 0.3, while for x > 0.4 it increases 
monotonically. Theoretical calculations in the framework 
of the RPA reproduce the experimental data reasonably 



well, both qualitatively and quantitatively. In contrast, 
the MFA not only overestimates the Curie temperature, 
but also exhibits a monotonic increase in the whole con- 
centration range, in contrast to the experiment. Accord- 
ing to Stearnsjii the larger T c of Cu2MnSn compared 
to that of Ni2MnSn is due to stronger localization of d- 
electrons, which enhances the Curie temperature. The 
flat minimum is explained by the competition of this 
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effect with the weakening of this interaction for larger 
interatomic distances (larger lattice constant) with in- 
creasing Cu-content.— It is obvious that this behavior is 
due to AFM-like exchange interactions of the 3rd shell. 
For low Cu-concentration this interaction somewhat re- 
duces the RPA sum, and for higher Cu-content its effect 
is compensated by increasing FM-like interactions. 

The calculated Curie temperature for Cu2MnSn is also 
shown in Fig. [7] Experimental values are not discussed in 
Ref§ There is some controversy in the literature: values 
530 K and 630 K can be found, although the lower one is 
sometimes questioned. The calculated RPA value 560 K 
lies in the range of the reported experimental values. 

The results for N^MnSn alloy can be compared with 
other theoretical work.— The experimental values of the 
Curie temperature range from 328 K to 360 K. The Curie 
temperature estimated in Refl39l using the MFA and the 
FM reference state is 323 K, if only the Mn-Mn exchange 
integrals are considered. Our related test value is 315 K. 
The Curie temperature rises to 360 K if Ni-Mn exchange 
interactions are also included^ Our mean-field value 
based on the DLM reference state is 373 K. We would 
like to point out that good agreement of the Curie tem- 
perature obtained in RefjH from the FM reference state 
with the experiment is probably fortuitous, as the RPA 
would lower this value. The present RPA value based 
on the DLM reference state (322 K) agrees well with the 
experiment. 



2. Extraction of exchange interactions from experiment 

Noda and Ishikawa^S have extracted up to 8 nearest- 
neighbor (NN) exchange interactions for Ni2MnSn and 
Pd2MnSn from the fit to measured inelastic neutron spin- 
wave scattering data. Taking care of the fact that the 
authors' exchange integrals do not include spin moments 
and different definitions of the prefactor in the Heisen- 
berg Hamiltonian, we compare in Fig. [8] our calculated 
results with the values extracted by Noda and Ishikawa.— 
The general agreement between the calculated results 
and the extracted values from the fit appears to be good. 
It should be noted that calculated exchange integrals are 
not limited by the distance 2 a (8 NN) as the fitted ones. 
In fact, saturation of results with respect to the number 
of shells is not achieved for 8 NN interactions (see Fig. [5]). 
It is also clear from Table I of Refl40i that the fitted val- 
ues can change significantly depending on the number of 
exchange integrals used for a fit. For example, in case 
of Ni2MnSn the 2nd neighbor interaction is larger than 
that of the 1st neighbor for the 8 NN fit, and smaller for 
the 6 NN fit. 

Ideally, one expects the fitted values to decrease in 
their sizes. This is not so for a relatively large J(8) in Ta- 
ble I for Ni2MnSn 40 which, for an optimal fit, should be 
of the smallest magnitude. This indicates that the fitted 
values should not be taken too literally. Seemingly per- 
fect agreement between the measured magnon spectrum 
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and that calculated from the fitted values for Pd2MnSn 
m RefS does not guarantee similar agreement in other 
related magnetic properties. 

We also estimate the spin-stiffness constants -D s tiff for 
Ni2MnSn and Pd2MnSn. The spin-stiffness is a property 
of the ferromgnetic state, i.e., the state at T=0 K. We 
have therefore estimated -D s tiff from calculated Mn-Mn 
exchange integrals in the FM state. The real-space ex- 
pression for Dstiff is non-convergent, and a special treat- 
ment is needed for its determination. 31 Estimated val- 
ues of L> s tiff for Ni 2 MnSn and Pd 2 MnSn are 160 ±25 
(150 ±10) meV.A 2 and 130 ±25 (98 ±10) meV.A 2 , 
respectively. The values in brackets are experimental 
values^ and the agreement between the theory and ex- 
periment is acceptable. The theoretical error bars were 
estimated from differences in calculated -Dstiff for various 
sets of damping constants^ 

Finally, corresponding results for (Nii_ x ,Pd K )2MnSn 
alloy are shown in Fig. [9] The experimental values^ , of 
the Curie temperature lie in the range 190-340 K. In line 
with the concentration dependence of the exchange in- 
tegrals (see Fig. [5J one expects a linear decrease of the 
Curie temperature with the Pd-content. The calculated 
RPA results show an almost linear decrease of the Curie 
temperature found in the experiment, although the de- 
crease is smaller than what is observed experimentally. 10 
The MFA values are systematically above the RPA as 
well as the experimental results, but do show the lin- 
ear decrease of the Curie temperature. Some disorder 
in Heusler alloys, even for stoichiometric compound, is 



rather frequent. In the case of Pd 2 MnSn alloy one can 
speculate about some disorder between Mn- and Sn-sites, 
which could reduce the calculated Curie temperature 
(due to antiparallel alignment of Mn-moments) . Such a 
study would require a generalization of the RPA to ran- 
dom systems, which is beyond the scope of the present 
article. 



D. Resistivity 

In general, there is very little known about the resis- 
tivity of Heusler alloys, either experimentally or theoret- 
ically. In magnetic alloys there are three contributions to 
the resistivity: a temperature-independent part, residual 
resistivity due to the alloy disorder and other defects, 
and two temperature-dependent terms due to electron- 
phonon and electron-magnon scatterings. 



1. Residual resistivity 

The residual resistivities are calculated using FM refer- 
ence state, appropriate for zero temperature. The calcu- 
lated resistivities for (Nii-^T^MnSn, T=Cu, Pd are 
shown in FigfTUJ The result, an almost parabolic con- 
centration dependence of residual resistivities (Nordheim 
rule)^ for both alloys, is in agreement with our discus- 
sion of the character of the disorder at the Fermi en- 
ergy (see Sec. IIII A"|) . namely its weakness. Much smaller 
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resistivity of (Ni,Pd)2MnSn alloy is due to the site off- 
diagonal disorder, whose effect on the resistivity is weaker 
than that of the diagonal disorder which dominates in 
case of (Ni,Cu) 2 MnSn. 



2. Temperature- dependent resistivity 

The contribution due to phonons is relatively well 
understood^ and it is small in Heusler alloys for temper- 
atures around the room temperature! 21 ' 22 The most im- 
portant contribution to the temperature-dependent resis- 
tivity is due to magnon scattering or the spin-disorder re- 
sistivity. A simple theoretical model of the spin-disorder 
induced resistivity of magnetic alloys was developed some 
time ago by Kasuya42- This theory explains the experi- 
mental facts that for temperatures higher than T c the 
resistivity is almost constant, while for temperatures be- 
low T c the resistivity decreases and it is equal to the 
residual resistivity at zero temperature. The increase of 
the resistivity with temperature (above T — and below 
T c ) is due to the increasing amount of the spin disor- 
der. The spin-disorder itself is best characterized by the 
spin-spin correlation function, 43 which can be determined 
from first-principles using known exchange integrals (see 
e.g. Ref|4J). Here we adopt a less ambitious approach. 
For temperatures above T Cl we assume that a reasonable 
model of spin-disorder is simply the DLM model, which 
assumes that the local magnetic moments on atoms are 
oriented randomly with equal probabilities in all direc- 
tions. The net moment is zero and the spins are un- 
corrected (the spin-spin correlation function is zero in 
this case). The DLM model is formally equivalent to a 
random alloy problem^ and the TB-LMTO-CPA code 
employing the Kubo Greenwood formula can be used con- 
veniently to compute the resistivity. 

Schreiner et al£^ studied the temperature depen- 
dence of resistivity of a series of Heusler alloys including 
Ni 2 MnSn and PdaMnSn^l The resistivities calculated for 
the DLM reference states for (Ni,T) 2 MnSn (T=Cu,Pd) 
alloys are presented in Fig. QT] We remind the reader that 
this resistivity accounts only for the spin-disorder and 
does not include the part due to scattering from phonons. 
We observe an increase of the resistivity due to spin- 
scattering for Cu- impurities in the host Ni 2 MnSn and a 
decrease for Pd impurities. Fig. ITT1 reflects the decreasing 
degree of localization (correspondingly, increasing Mn-X 
hybridization) along the sequence Cu-Ni-Pd, which leads 
to a decrease of the relative strength of the magnetic dis- 
order, and therefore resistivity, along the same sequence. 
The corresponding total resistivities measured in the ex- 
periment for the ordered Ni 2 MnSn and Pd 2 MnSn alloys 
are about 75 /lOhm.cm and 50 /lOhm.cm, respectively. 
The experimental values obtained after subtraction of the 
contribution from phonon scatterings are 47 pOhm.cm 
and 22—30 ^tOhm.cm. These values compare reason- 
ably well with the values 50 //Ohm.cm and 37 /lOhm.cm 
obtained from the present simple theory. In particular, 
larger value of the spin-disorder induced part of the resis- 
tivity in Ni 2 MnSn as compared to Pd 2 MnSn is correctly 
reproduced. 

We present a simplified treatment of the temperature 
dependence of resistivity due to spin disorder for tem- 
peratures below T e j 42 ' 43 We use two simple approaches, 
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based on (i) simulation of the temperature dependence 
through the uncompensated DLM (uDLM) modelr&^ 
and (ii) assuming experimentally observed quadratic 
temperature dependence of the resistivity due to spin- 
disorder in combination with calculated values of T c and 
the resistivity above the Curie temperature. 

The results of the first model are shown in Fig. [T^l For 
a given temperature below T c we identify the magnetiza- 
tion of the system with that corresponding to the uDLM 
model with a specific ratio of the parallel and antiparallel 
spins. We calculate the resistivity for this ratio using the 
conventional Kubo-Greenwood approach. If there are no 
antiparallel spins in the the system we have the FM state 
corresponding to the temperature T=0 K. On the other 
hand, for an equal number of parallel and antiparallel 
spins we have the DLM state with zero total magnetiza- 
tion, corresponding to resistivity at and above the Curie 
temperature. 

Optimally one can determine the magnetization of the 
system from the present Heisenberg Hamiltonian using 
statistical mechanical methods. Alternatively, one can 
adopt some models. We used two, namely the tem- 
perature dependence of the magnetization as obtained 
from the hyperfine field (HF) measurements on Sn atoms 
in Ni2MnSn^£ and that obtained from similar measure- 
ments on Fe-atoms in bcc-Fe4£ From the reduced mag- 
netization M/Mq at a given reduced temperature T/T c 
we determine the ratio of parallel and antiparallel spins 
for the uDLM model, and associate the corresponding 
reduced resistivity p/p c with the reduced temperature 
T/T c . For a typical ferromagnetic-f-i-paramagnetic tran- 
sition the reduced magnetization curve close to T/T c =l 
has a universal (system-independent) form. So our re- 
sults for p/ p c obtained this way should be more reliable 




T/T 



FIG. 12: (Color online) Reduced resistivity p(T)/p c vs re- 
duced temperature T/T c (the suffix c relates to the Curie 
temperature) for two models of the temperature dependence 
of the magnetization: (i) as the temperature dependence of 
the Sn hyperfine field in Ni2MnSn (squares, Ref \4w ■ and (ii) 
as the temperature dependence of the Fe hyperfine field in 
bcc-Fe (circles, Ref l47h . The experimental values including 
the phonon part (triangles) are from. Refl2ll. The dashed 
line is a simple quadratic curve, showing the magnon-part of 
the reduced resistivity. Being in the reduced form, it is ac- 
tually free of the constant B: p/p c = (T/T c ) 2 . The straight 
line is used as a guide for eyes. The inset shows the variation 
of the reduced magnetization M/Mq as a function of T/T c 
determined from Refl4r3 (squares) and Refl47l (circles). Mo is 
the saturation magnetization. 



in this range. A couple of comments are in order at this 
point. The HF field at Sn atoms in Ni2MnSn is not di- 
rectly proportional to the magnetization. Several models 
for the dependence of the HF field on the magnetiza- 
tion are discussed in the literature and one such model 
was used in Ref|46| by Gavriliuk et al. The seemingly 
linear decline of magnetization with temperature around 
T = contradicts the Bloch T 3 / 2 law, which a typical 
ferromagnetic material is expected to obey. It is in view 
of this and the uncertainties involved in relating the HF 
field at Sn atoms to the magnetization that we include 
the second model based on bcc Fe as the prototype of 
a ferromagneticoparamagnetic transition. The actual 
magnetization vs. temperature curve for Ni2MnSn is ex- 
pected to lie between the red (squares) and the green 
(circles) curves shown in the inset of Fig. [121 A typical 
convex form of the temperature dependence of the re- 
duced resistivity is well reproduced (see Fig. Q2]) for both 
cases. The experimental curve including the phonon con- 
tribution lies in between these two simple model curves. 
For comparison, in Fig. [T^] we also include the quadratic 
form pj p c — (T/T c ) 2 (dashed curve), based on the ex- 



12 




1 1 i i i i i i i i i i 

100 200 300 400 500 

r(K) 



FIG. 13: (Color online) The temperature dependence of the 
resistivity determined as described in the text for Ni2MnSn 
(triangles) and Pd2MnSn (circles). The empty and full sym- 
bols denote the calculated and experimental! results, respec- 
tively. 

perimental observation that the magnon resistivity has a 
simple quadratic dependence on temperature^ The dif- 
ference between the quadratic form (dashed curve) and 
the experimental curve is the contribution due to the 
phonons and we see that it is quite small. 

In Fig. Q2] we display theoretically calculated resistiv- 
ity assuming the form 21 p(T)= p Q + AT + BT 2 . The 
quadratic term is due to the spin disorder ; 42 ! 43 the linear 
one is the phonon contribution (which is valid with ex- 
ception of very low temperatures) j& and po is the value of 
the residual resistivity at T=0 K, which is due to all the 
defects present in the sample. While the values of po and 
the coefficient A were taken from the experiment^ the 
coefficient B was determined from the calculated T c and 
the resistivity in the DLM state. We obtain good agree- 
ment between the theory and experiment, in particular 
for Ni2MnSn. It should be noted that some offsets of 
the calculated curves as compared to the experiment are 
due to differences in calculated and measured T c which is 
larger for Pd2MnSn. Theoretical estimates of the coeffi- 
cient B are 4.79- and 5.2xlO"VO.cm K~ 2 for Ni 2 MnSn 
and Pd2MnSn, respectively. The corresponding experi- 
mental values are 3.94- and 6.1xlO"VO .cm K 2 . Con- 
sidering the fact that B represents the second derivative 
of resistivity with respect to temperature, this agreement 
is good, as evidenced by Fig[T3J 

IV. CONCLUSIONS 

We have studied the electronic, magnetic, thermody- 
namic, and transport properties of quaternary Heusler 



alloys (Ni,T) 2 MnSn (T=Cu, Pd) by means of first- 
principles density functional method. In agreement with 
experiments, magnetic moments per formula unit depend 
only weakly on the alloy composition for both alloys, hav- 
ing values around 4 /xb- Exchange interactions were de- 
termined using the DLM reference state, which assumes 
no a priori magnetic ordering in the system. The Curie 
temperatures were estimated using RPA applied to the 
non-random Mn-sublattice. The alloy disorder strongly 
influences values of exchange integrals in (Ni,Cu)2MnSn 
alloys, while only weak concentration dependence of ex- 
change integrals is found for (Ni,Pd)2MnSn alloys. Con- 
sequently, a linear decrease of the calculated Curie tem- 
perature is obtained in the latter alloys in agreement with 
the experiment. A more complex concentration depen- 
dence of exchange integrals in (Ni,Cu)2MnSn alloys, in 
particular different behaviors of Ni-rich and Cu-rich al- 
loys, can be ascribed to the 3rd NN exchange integrals. 
This result also confirms a model of Stearns^ on the 
relevance of the first three exchange integrals for mag- 
netic properties of Heusler alloys. In fact, a much larger 
number of exchange integrals is needed to obtain reason- 
ably stable results in ab initio theoretical studies. This 
is in striking contrast with the related semi-Heusler al- 
loy (Cu,Ni)MnSb,— where the exchange interactions are 
strongly damped due to their halfmetallicity. 

The residual resistivities obey the weak-scattering 
Nordheim rule. This is due to the fact that strong disor- 
der found e.g. in (Ni,Cu)2MnSn alloys influences energy 
states far from the Fermi energy relevant for electronic 
transport. Using a simple model for the spin-disorder, 
we have estimated the temperature dependent resistivity 
at temperatures above T c . Reasonably good agreement 
with experimental results is found for calculations which 
assume that the resistivity above T c is essentially cap- 
tured by the DLM model. Using this value of resistivity 
and calculated Curie temperatures, good agreement with 
the experiment was obtained also for the temperature de- 
pendence of the resistivity below the Curie temperature. 
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